Few-Body Systems 0, 1-11 (2008) p 

bystems 



© by Springer- Vcrlag 2008 
Printed in Austria 



^ QCD Phenomenology of Static Sources 

o ' 
o. 

fNl . Gunnar S. Bali 

> 

Department of Physics & Astronomy, University of Glasgow, Glasgow G12 8QQ, UK 

Abstract. We discuss the spectrum of open string and point particle excita- 
tions in QCD with various source representations. Some general relations are 
introduced and lattice results presented. In particular we discuss the short- 
q ! distance behaviour, relate this to perturbation theory expectations and com- 

CN ] ment on the matching between low energy matrix elements and high energy 

Wilson coefficients, within the framework of effective field theories. 

o 

r*| ■ 1 Introduction 



a: 



We discuss the excitation spectrum of QCD in its static limit. This is quite 



D . amusing, in the context of a Light Cone Workshop. However, some of the results 

are instructive and of a very general nature, in particular those concerning the 
spectrum of QCD and the question of power divergences and renormalons. For 
instance the Wilson-Schwinger line appears at prominent places. We shall see 
that in a scheme without a hard cut-off, such as dimensional regularisation (DR), 
such objects require special attention. Wilson lines do not only appear in the 
static limit, within the framework of heavy quark effective field theories, but for 
instance also within light cone parton distributions pQ, if one wishes to define 
them in a gauge invariant (and hopefully path-independent) way [2J G3 II]- 

In general, a Wilson-Schwinger line of length I in Euclidean space-time within 
a correlation function will result in a term oc exp(— Al), with some self-energy A. 
In DR any perturbative contribution to A vanishes. This is different in schemes 
with a hard momentum cut-off, such as lattice regularisation. In this case A 
contains a contribution that is proportional to the inverse lattice spacing a -1 
and which can be expanded in powers of a s . Within physical observables such 
power terms obviously have to cancel. 

To translate A from the lattice into the on-shell (OS) scheme the perturbative 
expansion of A has to be subtracted, replacing the a -1 term by a renormalon 
ambiguity. For physical observables to remain unaffected, the renormalon of A 
in the OS scheme has to cancel against a similar renormalon associated with 
a different term. One nice thing about a hard cut-off is that the structure of 
power terms within such a regularisation scheme is indicative of the renormalon 
structure that one will encounter in a DR calculation. One way of elucidating 
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this connection has been explained above: the price of cancelling a power term 
is a renormalon. Vice versa it is possible to remove (part of) such a renormalon, 
by introducing a scale dependence and replacing it by a power term. One such 
continuum scheme, the RS scheme, has been suggested in refs. El - 

Here I discuss the spectrum of QCD in the static limit. The simplest case 
would be that of a static-light meson. In the continuum limit the mass A of such 
a particle will diverge. Only level splittings are well defined. One can however, 
relate this mass to the mass of the physical B meson, within Heavy Quark 
Effective Theory (HQET). In this case the power term/renormalon of A will be 
compensated for by a similar contribution to m&. 

This article is organised as follows: in Section |2] we introduce the relevant 
correlation functions and discuss some fundamental relations. Then in Section 01 
we survey some aspects of the large distance behaviour and string breaking, 
before we discuss short distance relations in Section |1J 



2 Static QCD strings and particles: the broad picture 

We discuss the spectrum of QCD with one or two static charges, where we 
restrict ourselves to charges in the fundamental representation, 3, and the adjoint 
representation, 8. The former can be motivated as representing the starting point 
for realising a heavy quark within an effective field theory such as HQET, non- 
relativistic QCD (NRQCD) or potential NRQCD (pNRQCD) 00. The adjoint 
representation can be used to describe a gluino, within a heavy-gluino expansion 
(should a gluino exist and be sufficiently stable). 
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Figure 1. Correlation functions in the static sector of QCD. 



We schematically depict the correlation functions C(r,t) that we investigate 
in Fig. ^ The states are created at time and destroyed at time t. Straight 
spatial lines are parallel transporters, straight temporal lines static propagators 
(Wilson-Schwinger lines) in fundamental or adjoint representation. Curved lines 
represent sea quarks and gluons. The corresponding (in most cases r-dependent) 
energy levels can be extracted in the limit of large Euclidean times, 

T-i/ \ r dhxC(r,t) , , 

E(r) = - km . (1) 

t—>co at 
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We denote the energy level of a static-light meson by, 

M B = m b + A + 0(l/m b ), (2) 

where both the quark mass m b and the binding energy (or static energy) A 
are scheme- and scale-dependent, in reflection of the fact that asymptotic states 
cannot carry strong charge. Mb however is a physical observable. Note that 
m b does not appear as a scale in the static limit. I introduce it nonetheless for 
pedagogical reasons. Taking account of light quark spin and angular momentum 
the spectrum of A is arranged in spin-doublets (n + where for instance ^ 
corresponds to / = and ^ + to / = 1. The presence of such spin doublets has 
found an exciting application, with the recent discovery of two new D s states (HI 
E3 ini E]- Once l/m b corrections are added, the heavy quark spin has to be 
taken into account and the \ ground state will further split into a pseudoscalar 
B meson and a vector B* meson. 

The next step is an open string state, between a fundamental source and 
anti-source separated by a distance r. In this case the inclusion of l/m b cor- 
rections is essential, even at leading order as without a kinetic term all sorts of 
pathologies can arise which are related to the fact that heavy quarks at relative 
speed decouple from each other. The more conceptional results presented here 
are not affected by this problem. One encounters a ground state E s {r) as well 
as hybrid excitations En(r), due to radial and/or orbital gluonic excitations. 
Obviously, sea quarks screen these potentials and, 

lim E s (r) = 2Mb, (3) 

r — >oo 

which is known as "string breaking" . Note that the heavy quark mass mj cancels 
from the combination, E s {r) — 2Mb- A similar relation will apply to Eh, which 
will also decay, depending on the quantum numbers H, into a pair of B mesons. 

Within the framework of pNRQCD which amounts to a double expansion 
both in the relative quark velocity (NRQCD) and in the distance r, one can 
classify the short distance behaviour. In the static limit the only scale, apart 
from a typical non-perturbative QCD scale Aqcd, is the distance r and one can 
identify [Jj], 

E s (r) = 2m b + V s {r) + 0{r 2 ) = 2m fc L + £ s L (r), (4) 
E H {r) = 2m b + V {r)+A H + 0(r 2 )=2m b L + E H (r), (5) 

where both, singlet and octet potentials V s (r) and V (r) are calculable in pertur- 
bation theory for distances, r <C Aq 1 cd . Ah stands for the mass of a gluelump 
(see below). Again, m b , V s (r), V {r) and Ah are scheme and scale dependent. 
The new ingredient now is that the definition of the infrared gluelump energy Ah 
carries an ultraviolet renormalon ambiguity that has to be compensated for by 
V Q . The ultraviolet ambiguity of V Q in turn is the same as that of V, and related 
to the definition of the quark mass m b . On the other hand one can also calculate 
these energy levels non-perturbatively in lattice simulations and identify these as 
the ground state static potential, E^(r) (which is in the s = 17+ representation 
of the cylindrical group -Dqo/i) and hybrid excitations, E H (r). 
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The second equalities within Eqs. (pjj) and © above are only correct up to 
lattice artefacts which for the lattice actions used here will be proportional to 
a 2 AQ CD and a 2 /r 2 . The a 2 /r 2 correction means that these relations do not hold 
in the limit r — > 0. In particular, unlike the continuum E s , E~*(0) = does not 
diverge and E H (0) = A H . Note that in lattice schemes both m& and V s contain 
power divergences oc a" 1 which cancel each other. As r < a the above continuum 
limit interpretation is not valid anymore and E s (0) = V^z,(0) = 0. The lattice 
"1/r" term within V S) l assumes a finite value, oc a -1 , at r = 0, that exactly 
cancels the power term of the static energy. 

If we subtract the above equations from each other we obtain a consistent 
picture, 

E H (r) - E s {r) = V s {r) - V a (r) + A H + 0(r 2 ) = E H {r) - E^{r). (6) 

The second equality only applies modulo lattice artefacts, i.e. for r 3> a. This 
combination does not contain the leading renormalon and is also power term 
free. 

As indicated in Fig. ^ in the limit r — > the non-perturbative energy levels 
calculated on the lattice will either correspond to gluelumps or to glueballs, at 
least if we neglect sea quarks for the moment. In fact, as gluelumps are accom- 
panied by power divergencies, at sufficiently small lattice spacing, gluelumps will 
always form the ground state within each J PC sector. Starting from these glue- 
balls and pulling the static sources apart (1 © 8 = 3 (g> 3*) these states can be 
identified as the ground state in the fundamental open string sector plus glueball 
scattering states. If we start from a gluelump, increasing the distance will lead 
us to hybrid energies (or hybrid energies plus glueballs). At r > the symmetry 
is reduced from 0(3) <8> C to AxA and each gluelump state will in fact be ap- 
proached by more than one hybrid energy level. Note that gluelump masses are 
scale and scheme dependent while gluelump mass splittings are universal. 

We have restricted the discussion of gluelumps and hybrid potentials to the 
case without sea quarks. Note however that the presence of (massive) sea quarks 
will not change the situation conceptionally. We have already discussed the break- 
ing of the ground state string above and the same will happen in the hybrid 
sector. Sea quarks will result in BB potentials, in addition to the conventional 
hybrid potentials. Their presence will provide us with an alternative possibility 
of screening of a static octet charge by quark and antiquark and this will give 
rise to additional scattering states. The level orderings will somewhat change 
and the excitation spectrum will become more dense. However, the pNRQCD 
multipole expansion also provides the framework for classifying this situation 
and moreover the renormalon and power term structure will remain unaffected. 

Finally, the gluelump energy Ah is an object very similar to the binding 
energy A of a B meson. Imagine a heavy gluino with mass rrig. This will be 
screened by the gluons within the QCD vacuum and only visible as a bound 
state glueballino with mass, 

M d = m~ g + A B + 0(l/m~ g ). (7) 

The lightest glueballino will be the magnetic one, H = l" 1 = B. The scheme 
and scale dependence of Ab is required to cancel that of m§. In a similar way 
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in which heavy- light mesons are related to quarkonia at r — » oo, glueballinos are 
related to gluinonia, bound states of two heavy gluinos: 

E A (r) = 2m~ g + V A , s (r) + 0(r 2 ) ™ 2M^, (8) 

where VA,s( r ) is the singlet potential between two adjoint sources. As 8 © 8 = 
Iffi8©8ffil0ffi 10* © 27, at r — > contact can be made between static hybrid 
gluinonium energy levels and gluelumps, this time not only in representation 8 
but also in 10 and 27: a whole tower of states exists, and each new representation 
introduces a new renormalon/power term. 

Unfortunately, no numerical results on hybrid excitations of non-fundamental 
QCD strings or of higher representation gluelumps exist thus far. However, sin- 
glet open string excitations in representations even larger than the adjoint one 
have for instance been calculated in ref. ^3] an d the closed string spectrum in 
ref. [Hj. 

3 String Breaking and nib 

Attempts to resolve string breaking in lattice studies of QCD with sea quarks 
have a long history. While everyone knows that this effect exists, it is still a crucial 
benchmark for the capability of lattice calculations. Moreover, the dynamics of 
the mixing between open string and broken string states, which are the starting 
points for a description of quarkonia and heavy-light mesons, respectively, is non- 
trivial. This will elucidate strong quarkonium decay rates and other properties 
near threshold. 




2 4 6 8 10 12 14 16 18 2 4 6 8 10 12 14 16 18 

r/a r/a 



Figure 2. The lowest two energy levels in the QQ system with sea quarks. On the right hand 
side the BB content of the ground state is shown, in terms of a mixing angle. The results apply 
to n/ —2, m q ~ m s and a ~ 0.085 fm. 

First studies in toy models were performed as soon as in 1988 ^H]- Only some 
five years ago quantitatively satisfying results were obtained, first in 5t7(2)-Higgs 
models ^3 El and then for the adjoint representation string in SU(2) gauge 
theory ^3^2- The latter case corresponds to the decay of gluinonium into two 
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gluinoballs, discussed above. Several attempts on the string breaking problem in 
QCD with sea quarks were made j2UJ |^ but only very recently reliable 

results were obtained 123] • These are displayed in Fig. [21 String breaking takes 
place at around r c ~ 1.27 fm for light quark masses similar to that of the strange 
quark. An extrapolation to physical quark masses yields r c ~ 1.16 fm. The gap 
between the two states in the string breaking region is about AE « 50 MeV 
and we are able to resolve this with a resolution of 10 standard deviations! One 
might expect lighter sea quarks to extend the pion exchange related "bump" in 
the excited energy level and mixing angle towards larger distances and also to 
further broaden the gap AE. 

2.5 
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Figure 3. The binding energy A from the QQ potential (solid diamonds) and from static-light 
mesons (open diamonds) as well as the perturbative conversion into the RS scheme at different 
orders. 

We have discussed that E s (r) — > 2Mb for r — > oo. This is not the case in the 
quenched approximation, however, still at some distance r c = 1.13(7) fm [6j 121 j 
we will find ^E s (r c ) = Mb- On the lattice we have, 

\E^(r c ) = \E s (r c ) - m biL = 3p ot . (9) 

One can determine both, the binding energy from the static potential, Ap Qt (a) = 

^Eg(0.5 fm; a) + A, up to a constant A and A L (a) directly from the static- 
light system [251 1261 1271 112j . The latter values are less accurate since a chiral 
extrapolation in the light quark mass is required. There are also inconsistencies 
between different A data sets obtained by different groups at the coarser lattice 
spacings. 

The leading overall renormalon ambiguity cancels from the running of A L 
from one scale vj = a -1 to another. In Fig. EJwe plot the results (diamonds), 
together with the NNNLO perturbative expectation, expanded in terms of 
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a s (3.9GeV). We find excellent agreement between perturbation theory and the 
lattice data, down to energies as low as Vf < 1 GeV. We also translate the results 
into the RS scheme jH]. At large scales (small lattice spacings) the power term 
becomes large and hence a conversion between schemes requires very accurate 
perturbative coefficients. At small scales this requirement is less demanding but 
perturbation theory obviously becomes less convergent. AT NNNLO the optimal 
accuracy can be obtained around Uf « 4 GeV. 

By subtracting the power term from A one obtains the binding energy in the 
OS scheme. This is scale independent but contains a renormalon ambiguity. If we 
are interested in extracting the b quark mass in the MS scheme the renormalon 
above will cancel against the one that arises from converting the OS b quark 
mass m,b t os{ a ) = Mb — A (a) into the MS mass [2H1- We use the RS scheme 
as an intermediate scheme in this conversion and obtain [H], 

m bMs( m b,Ms) = t 4191 ± 29 ( stat -) ± 47 ( th ') ± ^MS)] MeV - ( 10 ) 
The theoretical error includes l/nib corrections. In addition to the errors dis- 
played we expect an 0(100 MeV) quenching uncertainty. 



t 1 1 r 




j i i L 



0.5 1 1.5 2 2.5 

r/r 

Figure 4. Hybrid energies 29 , in comparison with the gluelump spectrum, extrapolated to the 
continuum limit 30. (circles, left-most data points). The gluelump spectrum has been shifted by 
an arbitrary constant to adjust the l" 1 state with the II U and E~ potentials at short distance. 
The lines are drawn to guide the eye. 

4 Gluelumps and Hybrids 

At short distances one would expect more than one A° v hybrid to approach the 
same J PC gluelump energy level where J > A and rj = PC since -DqcA C 0(3)<8>C 
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Qualitatively this can be verified in Fig. |1] where we compare hybrid energies 
obtained in ref. j^H] with the gluelump spectrum of ref. |3U| . The dashed lines 
correspond to our expectations. Only the S^" level cannot be disentangled from 
a glueball scattering state. Everything is displayed in units of tq ~ 0.5 fm. 




0.25 0.5 0.75 1 1.25 1.5 



r/r 

Figure 5. Splitting between the IJ U and the energies. (V ,rs — V s ,rs)(»") +/1§ S is plotted at 
tree level (dashed line), one- loop (dashed-dotted line), two loops (dotted line) and three loops 
(estimate) plus the leading single ultrasoft log (solid line). 

We calculate the lowest two hybrid potentials E H (r) with H = IJ U and H = 
U~ which both will approach the lightest (magnetic) 1 h gluelump. According 
to Eq. ©, the hybrid energies are the non-perturbative generalisation of the 
octet potential of the perturbative multipole expansion. pNRQCD predicts the 
leading order difference between these two levels to be oc r 2 , which we are able 
to verify. 

There are now two strategies of determining a gluelump mass A B - Either one 
can directly calculate it at a given lattice spacing a on the lattice and subse- 
quently convert it into another scheme or one can compute the power term free 
combination E n (r) — E^(r) at r > a and extrapolate the result to the con- 
tinuum limit. Subsequently, one can then follow Eq. JHJ), subtract V s (r) — V {r) 
and obtain A B . The result of a continuum limit extrapolation of the difference 
is displayed in Fig. together with the expectation Eq. © to different orders 
in perturbation theory, using the RS scheme jSHH]. The nj = QCD coupling 
has been taken from ref. jHJ and hence the only free parameter in the fit is the 
gluelump energy, 

A B S {1 GeV) = [887 ± 39(latt.) ± 83(th.) ± 32(%g)] MeV. (11) 

We obtain a compatible result, following the first strategy outlined above, 
perturbatively converting the lattice gluelump data of ref. |30| at finite lattice 
spacings, A B {a), into the RS scheme, in analogy to our determination of the 
static-light binding energy. The second lightest gluelump is the electric one (1 ) 
which is about 350-400 MeV heavier than the magnetic gluelump. 
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5 Conclusions 

Relations between static energy levels at short and large distances (string break- 
ing) have been reviewed. For the static energies NNLO/NNNLO MS perturba- 
tion theory works very well down to energies of less than 1 GeV, once the leading 
renormalon has been accounted for. Static Wilson-Schwinger lines in the funda- 
mental and adjoint representations give rise to masses (the binding energies A 
and gluelumps, respectively) that are scale and scheme dependent. This has im- 
plications with respect to QCD vacuum models and condensates. In particular, 
when combining perturbative Wilson coefficients with non-perturbative matrix 
elements these have to be defined in the same scheme and at the same scale to 
enable renormalon cancellation. 
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